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We study when a non-local unitary operation acting on 
two eí-level systems can probabilistically simulate another one 
when arbitrary local operations and classical communication 
are allowed. We provide necessary and sufíicient conditions 
for the simulation to be possible. Probabilistic interconverta- 
bility is used to define an equivalence relation between gates. 
We show that this relation induces a finite number of classes, 
that we identify. In the case of two qubits, two classes of 
non-local operations exist. We choose the CNOT and SWAP 
as representatives of these classes. We show how the CNOT 
[SWAP] can be deterministically converted into any opera- 
tion of its class. We also calculate the optimal probability of 
obtaining the CNOT [SWAP] from any operation of the cor- 
responding class and provide a protocol to achieve this task. 
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Much of the attention in Quantum Information Theory 
[QIT] was focused in recent years on obtaining a quali- 
tative and quantitative description of the entanglement 
properties of pure and mixed states. Apart from the im- 
portant link of this problem to the bàsic principies of QIT 
and Quantum Mechanics, a proper understanding of en- 
tanglement is also expected to lead to new applications 
in quantum communication and quantum computation. 
As entanglement — mainly in the form of a maximally 
entangled state — is a key ingredient for many applica- 
tions in QIT, it became a focus of interest to consider 
entanglement as a valuable resource. In this context, it 
is important to know how to manipulatc entangled states 
in an optimal way, because this determines which QIT 
tasks can be done (with optimal probability of success) 
using a given entangled state. 

Quitc rccently, it was realized that entanglement of 
states is not the full story — one may also consider the 
entanglement properties of quantum dynamics. This is 
motivated in part by the facts that after all we deal with 
interactions in experiments, and the interactions allow 
us to create entangled states. This makes a proper un- 
derstanding of qualitative and quantitative entanglement 
properties of quantum operations highly desirable. In 
fact, first steps in this direction have been reported re- 
cently [|-[|]. However, even in the case of bipartite uni- 
tary dynamics — which may in some sense be considered 
as the analog of bipartite pure states — , only very few is 
known up to now. Important issues, such as the structure 
of non-local unitary operations as well as the possible 



local interconvertability of non-local unitary operations 
— i.e. the (probabilistic) simulation of an unitary oper- 
ation using some other one — are completely unexplored 
yet. Similarly as for pure state conversion, this last is- 
sue provides information about which kind of dynamics 
may be performed (with an optimal probability of suc- 
cess) given some unitary operation and allows us to put 
some structure in the set of non-local operations. 

Apart from the theoretical interest of such qüestions, 
they may also be of practical relevance. In particular, in 
many cases it is important to know whether a given re- 
source (i.e. a certain unitary operation) suffices to imple- 
ment a relevant task in QIT. For example, one may con- 
sider the the process of entanglement purification Q , one 
of the bàsic primitives for long range quantum communi- 
cation using quantum repeaters 0. The known schemes 
for entanglement purification Iqjrequire the possibility to 
implement CNOT operations ||] between two partides. 
The approach sketched above allows to decide whether 
a given unitary operation — e.g. produced by a weak 
interaction — is already sufíicient to implement a CNOT 
operation and thus entanglement purification. As entan- 
glement purification itself is already a probabilistic pro- 
cess, one can consider probabilistic simulation of the op- 
erations. In this letter, we provide a complete solution of 
the problem stated above. From our results follows e.g. 
that any (arbitrarily weak) bipartite unitary operation 
already allows to implement entanglement purification. 

We consider two (spatially) separated parties and non- 
local unitary operations U, U acting on two d-level sys- 
tems, each of them hold by one of the parties. We 
say that an unitary operation U can simulate the ac- 
tion of another operation U on an arbitrary input state 
probabilistically (equivalently, U can be converted to U), 
U — ^ U, if there exists a protocol using (stochastic) lo- 
cal operations and classical communication [(S)LOCC] 
Q and a síngle application of U to achieve this task. 
The most general simulation protocol in this context con- 
sists of appending local auxiliary systems, a sequence of 
(S)LOCC, followed by a single application of U and an- 
other sequence of (S)LOCC. We have to know when the 
protocol succeeds. In this case — which occurs with some 
non-zero probability of success-, the total action of the 
protocol on an arbitrary input state p — after tracing out 
eventually involved auxiliary systems — is given by U pU^ . 
Note that we do not allow for manipulations during the 
evolution (as is e.g. done in Hamiltonian simulation p]). 
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Such a situation appears if the process of interaction is 
inaccessible for some reason, e.g. because the interaction 
is taking place at a very short timescale, or the unitary 
operation is given by a hxed device (black box), which 
should then be used for some other purpose. 

One may also consider simulation of unitàries undcr 
other restricted sets of local operations, such as local 
unitàries (LU) or local operations without classical com- 
munication. For each kind X of local operations, one 
obtains in this way an equivalence relation based on in- 
terconversion of U and Ü, U U, which is discussed 
in detail in ]ÏC| ], Throughout this letter, we will consider 
probabilistic interconversion and allow for arbitrary local 
resources Q, i.e. SLOCC. 

The results we obtain are as follows: 

(i) We provide a necessary and sufhcient condition 
for probabilistic simulation of unitary operations 
U U. This allows us to obtain a complete, hier- 
archic classification based on the equivalence rela- 
tion induced by probabilistic interconvertability of 
operations, i.e. U # U. The number of classes is 
finite and for a(F d (g>(F d system, i.e. U £ SU(d 2 ), 
bounded from above by d 2 . 

(ii) For qubits, i.e. d = 2, we show that only three 
classes exist, which can be represented by the iden- 
tity, CNOT and SWAP respectively §. 

(iii) Also for the case of qubits, we explore the infer- 
nal structure of the classes and show that given 
the CNOT [SWAP] operation, one can determin- 
istically simulate any operation of its class. We 
also calculate the optimal probability to obtain the 
CNOT [SWAP] operation from an arbitrary opera- 
tion of the corresponding class and provide a prac- 
tical protocol to achieve this task. 

In order to obtain a necessary and sufhcient condition 
for probabilistic gate simulation (i), we make use of the 
isomorphism between non-local physical operations and 
states [Q. When applied to unitary operations, one hnds 
that to each unitary operation U S SU(d 2 ) corresponds 

a pure state \^u) e(F d ®d d given by 

|*tr) = ÜAiBxWlAiA* ® I$)biB 2 > (1) 

where \$)a x A 2 = Vv^S/Li I^aJ^)^ is a (local) max- 
imally entangled state. The isomorphism has a very sim- 
ple interpretation Q: (a) On the one hand, Eq. ([]]) telis 
us that one can obtain the state \^u) starting from a 
product state (systems A - B) given a single application 
of U. (b) On the other hand — as shown in jïj — given 
the state | W [7) , one can probabilistically implement the 
unitary operation U on an arbitrary input state p by per- 
forming suitable local measurements. We denote by 
the Schmidt number of the state Y&u), i-e. the number of 



non-zero Schmidt coefhcients. Recali that (c) a bipartite 
pure state \tp) can be transformed into another pure state 
\(j>) with non zero probability of success using SLOCC iff 
fhj, > (see e.g. Q). 

We can now state the following necessary and sufhcient 
condition for probabilistic gate simulation: 

U Ü iff n* u > n q , ü . (2) 

Using a sequence of local operations and the properties 
(a-c), the proof of Eq. (|J) can be summarized as follows: 
Necessity follows from the existence of the (probabilistic) 

local process \^u) U — > Ü ^> \^fj), which according 
to (c) implies that n^ u > n^ ü . Regarding sufhciency, we 

have: U H \^u) Q |%) ^> Ü. Note that the last rela- 
tion provides a protocol to achieve the simulation U U. 

It immediately follows that U can be interconverted 
into U probabilistically, U ^ U, iff nqi u = n>p & . We have 
that two operations U, U belong to the same equivalence 
class induced by this equivalence relation iff their corre- 
sponding pure states [^t/)!; \^ü) have the same Schmidt 
number. Since n^, u < d 2 , we have that at most d in- 
equivalent classes exist if U € SU(d 2 ). The classifica- 
tion is hierarchic, as unitary operations corresponding to 
states with a higher Schmidt number can simulate op- 
erations with corresponding states with a lower Schmidt 
number. 

In order to make this more explicit, we will now turn 
into the probabilistic conversion of two-qubit gates, i.e. 
d = 2 (ii). Wc cxplicitly obtain the corresponding classes, 
which turn out to be only three. As shown by Kraus et. 
al. H (see also B]), any two-qubit unitary operation 
U € SU (4) can be uniquely jï3] written in the following 
Standard form 

U AB = V A ® W B e- iH V A ® W B , (3a) 

3 

H = J2 H h H i = ® (3b) 

i=l 

n/4 > fií >H2 > |M3Í >0, (3c) 

where V A ,W Bl V A ,V B are local unitary operations. 
That is, up to local unitàries, any two-qubit uni- 
tary operation is given by the normal form e~ lH , 
which contains all the non-local content of the oper- 
ation. We denote by {|$í)}i=o,i,2,3 a basis of maxi- 
mally entangled states, where = iTj ® 11$) and 

|$) = 1/V5(|00) + 1 1 1>) . We have that \^ v ) = 
V A , ® W Bl ® Vl 2 ® Wl 2 e- lH ^\^) Al ^) Bl2 = LU 

Y,l=oai\®i)A 1 , 2 \®i)B lt2 - We used that l Al ® Va 2 \<S>) = 
Vj (8> 1a 2 |$) i n the first equality, while the second equal- 
ity is understood up to local unitàries and corresponds 
— up to some irrelevant phase factors — to the Schmidt 
decomposition of \^u) with corresponding Schmidt coef- 
ficients \oí\. We also used that e~ lH = J2 k=0 d%(rf ® of . 
The coefficients a% are given by 
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and we introduced the shorthand notation c u 



(4a) 
(4b) 
(4c) 
(4d) 

COs(/Xi), 

s Mi = sin(/ii). It is now straightforward to check that 
the Schmidt number is either 1,2 or 4. Thus three 
classes of two-qubit unitary operations under probabilis- 
tic local interconversion, U # U, exist: 
Class 1: riï&v = 1 : These are local unitary operations, 
with [i\ = [j>2 = A*3 = in Eq. (|^). One can choose the 
identity as a representative of this class. 
Class 2: n^s u = 2 : These are non-local unitary op- 
erations with fi\ 7^ and /i2 = (13 = in Eq. (||). 
The CNOT operation appears as a natural representa- 
tive of this class, as |*c/ CNOT ) = 75 (I 00 )a 1 , 2 |í > o)s 1 , 2 + 
|11)ai 2 2 ) i s a niaximally entangled state with 
Schmidt number 2. Note that the CNOT is up to lo- 
cal unitàries equivalent to an operation of the form (j^) 
with /ii = n/<i,(i2 — 113 = 0. 

Class 3: n^ v = 4 : These are non-local unitary opera- 
tions with /ii,/X2 7^ and arbitrary fj, 3 . The SWAP oper- 
ation appears as a natural representative of this class, as 
I*í7swap> = \®o)a 1 b 2 \$o)a 2 b 1 is a maximally entangled 
state with Schmidt number 4. Note that the SWAP is up 
to local unitàries equivalent to an operation of the form 
@ with ni = fi 2 = M3 = tt/4- 

Recali that any operations of class 3 can simulate op- 
erations of class 2 probabilistically, however the reverse 
process is not possible. This implies on the one hand that 
any non-local unitary operation can be used to simulate 
a CNOT operation probabilistically (and thus to implc- 
ment entanglement purification) , while the CNOT opera- 
tion can e.g. not be used to simulate e -1 *^ 1 ® ai +<T2 ®°' 2 ) 
with non-zero probability of success even for ícl. 

We now explore the internal structure of the different 
classes (iii). We first show that the representative of each 
class can be used to simulate any other operation within 
the class with unit probability. For the SWAP operation, 
this is trivial to see: Since the SWAP operation can be 
used to create two ebits of entanglement (see Eq. (Q)), 
one can use the first ebit to teleport a qubit from B to A, 
implement the desired operation locally in A and teleport 
the qubit back to B. Note that there is a whole set of 
operations within class 3 which can be deterministically 
interconverted into the SWAP operation. In particular, 
all operations of the form (|^) with /ii = /i 2 = tt/4 and 
arbitrary /X3 can create two ebits of entanglement and can 
thus similarly be used to deterministically implement any 
other operation. Note that these are the only two-qubit 
operations with this property. 

We now show that the CNOT can simulate any oper- 
ation of the form U(a) = e - * 00 ^® ^ with < a < tt/4 
deterministically. This is sufficicnt to obtain that any 



operation in class 2 can be simulated by the CNOT, as 
any operation in class 2 is equivalent to U (a) up to LU. 
Given an arbitrary input state \4>)ab = A\j) b , 

we use the following protocol: (a) We apply CNOT^^, 
where B is an auxiliary system at site B prepared in 
state |0). (b) We apply locally the operation U{a) BB . 
(c) We measure the particle B in the x-basis, i.e. we 
measure projectors corresponding to the states |±) = 
1/V2(|0) ± |1)). If we obtain the result correspond- 
ing to |— ), we apply 03 in A and B, otherwise we 
do nothing (the last step requires classical communi- 
cation). It is straightforward to check that the pro- 
tocol performs the desired gate on an arbitrary input 

Y,i,jhj\ï)A\j)B\Q)B ^ Y,i,jbij\i)A\j)B\ï)B ^ 
Ma(-l) ííl+Í41 )&ii|Í>A|j>J ' 



state: 



)B\l 



ls Í a (-l)^+^)b l] \ ï } A \j) B \+) B = U AB \4>)ab\+)b; 

In the remainder of this paper, we will consider op- 
timal simulation of CNOT [SWAP] given an arbitrary 
operation within the corresponding class, that is we are 
looking for the simulation protocol with highest possi- 
ble probability of success. We show the following: (a) 
If a protocol can probabilistically simulate the action of 
a unitary operation on all input states, then the prob- 
ability of success is independent of the input state. (b) 
We derive an upper bound for the probability of success 
for a particular input state, which — according to (a) — 
provides also a bound for the optimal protocol, (c) We 
provide a protocol which reaches the upper bound. 

Regarding (a), we use that any simulation protocol can 
be considered as a two-branch protocol with two possi- 
ble outeomes: The successful branch described by the 
operator M, where with some probability of success U is 
performed and the second branch where the simulation of 
U failed. We consider the successful branch and assume 
that the probability of success, p^ — tr(M|?/>)(?/;|M'ï) 
of the protocol depends on the input state \ip), where 

ÉtíOiily}- We have that M IV-> = y/P+UW)- us- 

r pïjU\ij), we have on the one hand 



ing that M\ij) ■■ 



Ei 



PipUj2ij a ij Ni) an d on the other hand 



M J2ij a ij\ij) = Htj a ij-JPïïU\ij). This implies that 
y/jty = ^Jpïj \fij and thus p^ = p is constant. 

(b) It follows from (a) that we can consider an arbitrary 
input state to derive an upper bound for the probability 
of success of the optimal protocol. As CNOT [SWAP] 
can create one [two] ebit(s) of entanglement when ap- 
plied to product input states (see Eq. (Q)), the opti- 
mal simulation protocol — when acting on a product in- 
put state — should also be capable of creating one [two] 
ebit(s) of entanglement. The first sequence of SLOCC 
does not change the product structure of the input state. 
Thus we have to consider the optimal probability to ob- 
tain a maximally entangled state with one [two] ebit(s) 
of entanglement using a single application of U acting 
on an arbitrary product input state followed by SLOCC. 
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Let |</>)a|x)b be an arbitrary product state, already in- 
cluding local auxiliary partides. We denote \4>í)a 1 a 2 = 
af 1 ® 1a 2 \(/>)a 1 a 2 and \xí)b x b 2 = crf 1 <g> ls 2 |x)s 1 s 2 , 
which are normalized, but not necessarily orthogonal 
states. Using the notation introduced between Eq. (|^) 
and Eq. (||) and assuming without loss of generality 
that U — e~ lH (the local unitary operations of a gen- 
eral U of Eq. (|J) can be absorbed into the definition of 
|<^>), |x) and the SLOCC performed afterwards), we ob- 

tain \xu) = U Ai b 1 \4>)a\x)b = J2k=i a>k\<f>k) A\Xk) b ■ Let 
bi be the ordered Schmidt coefhcients of \xu), bo > b\ > 
b 2 > b 3 . Using the results of (Ï^JÏ^l on optimal conver- 
sion of pure states, we have that the optimal probability 
to obtain a maximally entangled state of one [two] ebit(s) 
out of \xu) is given by p x = 2(6? + b\ + bf)[p 2 = 46§]. 
Note that for operations in class 2, the maximal prob- 
ability pi reduces to p\ = 2b\, since b 2 = 63 = in 
this case. Let \4>j)A be a normalized state for which 
K$;|0í}| = f° r i ¥"3- For \ t )a being an arbitrary state 
and pa = ^b 1 b 2 (\xu)(xu\) the reduced density matrix of 
system A, we have that ò§ < (ta\pa\ta)- For |ta) = \<t>3), 
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< (hÏPAlh) < |a 3 r|(03|03)| 2 < |a 3 | 



we obtain ò| 

This implies that p 2 < 4|a3 1 2 . A similar argument can 
be used to obtain a bound for p\ for unitary operations 
within class 2, however one has to take into account that 

which implies 



and thus obtains b\ < 



2s 



b 2 = 63 

Vi <2| ai | 2 - _.,„. 

(c) The optimal protocol is given by the one sketched 
in the proof of Eq. (§), i.e. U — > [*(/) -> 
\^u CNOTISW ap]) -> U C not[Uswap], where in the sec- 
ond step the optimal protocol for pure state conversion 
p2| , p^| is used. Note that the first and third step can be 
performed with unit probability of success, where the last 
relation follows from the fact that one [two] ebit(s) of en- 
tanglement are sufficient to implement a CNOT [SWAP] 
deterministically . We have that the conversion prob- 
ability to obtain a CNOT [SWAP] given an operation in 
class 2 [3] is Pl = 2\ ai \ 2 = 2s^ [p 2 = 4|a 3 | 2 ], which 
reaches the upper bound derived above and the protocol 
is thus optimal. 

In this letter, we derived necessary and sufficient condi- 
tions for probabilistic (inter) conversion of bipartite uni- 
tary operation and obtained a complete, hierarchic clas- 
sification. For two-qubit operations, we proved the exis- 
tence of three inequi valent classes, represented by Iab, 
CNOT and SWAP. We provided protocols to obtain the 
optimal conversion between the representative of each 
class and an arbitrary operation of the class. Note that 
the results presented in this letter can be extended to 
multipartite systems and one may also consider inter- 
conversion of operations under restricted classes of local 
operations |ÏCj . 
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